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We review the features of the bound state equation in large N scalar QCD in two dimensions, the 
't Hooft model, and compute the discrete hadron mass spectrum in this theory. We make the Ansatz 
that the scalar fields of this model represent spin zero diquarks and we estimate the minimum allowed 
mass for the first radial excitation of the lowest diquark-antidiquark scalar meson. The discussion 
is extended to the case of spin one diquarks. 



Introduction. In a recent paper on scalar meson dynamics it is shown how a satisfactory 
explanation of light scalar meson decays can be reached by assuming a dominant diquark-antidiquark 
structure, qq, for the lightest particles. The diquark q is a spin zero antitriplet color state. In 
first approximation the nonet formed by /o(980), ao(980), k(900), cr(500) is interpreted as the lowest 
qq multiplet. On the other hand, the decuplet of scalar mesons with masses above 1 GeV, formed 
by /o(1370), /o(1500), /o(1710), ao(1450), i4:o(1430), and likely containing the lowest glueball, is inter- 
preted in (Ij] as the lowest qq scalar multiplet (see also (3). The underlying hypothesis that the latter 
multiplet is not a radial excitation of the former has never been proven. 

In this note we provide evidence in support of this hypothesis by estimating the mass of the first 
radial excitation of the lowest sub-GeV qq scalar meson. 

We perform this calculation in large-TV QCD in (1 + l)-dimensions. This is a planar, linearly 
confining theory which admits a Bethe-Salpeter equation describing the discrete spectrum of qq bound 
states [3|. In this theory no orbital angular momentum excitations are possible since no rotation 
operator can be introduced: the discrete spectrum describes radial excitations. 

Our Ansatz is that the scalar fields of large N scalar chromodynamics in (l-l-l)-dimensions (SQCD2) 
can be thought as diquark fields. The corresponding Bethe-Salpeter equation of this theory should 
then yield the spectrum of tetraquarks qq particles. 

The bound state equations in spinor and scalar chromodynamics are respectively (see e.g. Q): 
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= ^^/i-) - I dyj-—^m (1) 

where the integrals are in the sense of the Cauchy principal value and 0(0) — 0(1) = 0. Here /i^^^ 
and /i^^-j are the mass squared eigenvalues of standard qq mesons, call them pions, and tetraquark 
qq mesons, sigmas. The parameters Mf and Ms represent the masses of the quark q and of the 
scalar quark q respectively. In this notation all the masses are adimensional parameters since we 
are rescaling them by the couplings g^: in two-dimensional scalar and spinor chromodynamics the 
couplings have dimension of mass. 

We cannot expect that the mass spectra in (s)QCD2 reproduce numerically the physical values of 
the masses of real pions and sigmas but we can assume that the regularities in the spectra of this 
kind of hadron-string models resemble those in the physical ones. To begin, we require that the ratio 
between the ground states /i((T)/M(7r) corresponds to the ratio between the lowest lying tetraquark 
object, the (t(500), and the lowest standard qq meson, the pion. We find that this is allowed by a 
family of values of the parameters {Mf, Ms). 

The first observed meson with the same quantum numbers of the pion is the 7r(1300) Q which is 
therefore a good candidate for the first radial pion excitation. Considering that all bound states are 
alternately even or odd under parity [6] and that the ground state for a qq particle has positive parity, 
the bound state corresponding to 7r(1300) is the second excited state in the spectrum /i(7r)- Then a 
simple calculation allows to estimate the mass of the first expected physical radial excitation of the 
cr(500) meson. 

It turns out that, spanning the set of allowed masses {Mf,Ms), the minimum value for the first radial 
excitation of the (t(500) tetraquark is at about 2600 MeV. This would predict an extremely broad 
state hardly identifiable experimentally. Moreover this result underscores that there is no possibility 
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that the above 1 GeV scalar mesons m the range 1350 1700 MeV could be excitations of the below 
1 GeV ones. This strengthens the hypothesis that in the latter mass region there is space only for the 
lowest qq multiplet (and likely a glueball). 

In the last section we extend the discussion to the case of heavy-light diquarks introduced to describe 
the X,Y,Z particles found by Belle and BaBar. In particular we will focus on the charged Z(4430) 
particle recently observed by Belle. As this is a spin one resonance, spin one diquarks have to be 
considered. We find that the bound state equation for scalar chromodynamics can be exploited also 
for a discussion of heavy-light axial diquarks. 

We start by illustrating the method used to solve Eqs. (|1I2[) and the numerical results obtained. 

Large N scalar QCD2- We will briefly review scalar chromodynamics in (1+1) dimensions and 
outline the derivation of Eq. A thorough discussion of Eq. ([T]) in spinor chromodynamics is found 
in Omitting color indices, the scalar chromodynamics Lagrangian density in light-cone coordinates 



X — 



V2 



(3) 



C = d^ip^d^^ - m\^^ + \{d^A+){d+A-) - ^A+(^t5_^ _ (a_^t)^) (4) 

^ V A 



reads: 



where the ^ index is /i = +, — , the gauge fields are in the adjoint representation of SU{N) whereas 
the scalar fields arc in the fundamental one and the light-cone gauge A'^ — A- — has been imposed. 
In this gauge no self-couplings of the gauge fields or seagull terms exist. The theory is asymptotically 
free and linearly confining. To make the self energy smooth at large A, we introduce I/Va" at each 
vertex. 

The emergence of a linearly confining potential can be observed by writing the equation of motion 
for the field A+: 

SiA, = -^^1- (5) 



which admits the following solution (use d-\x | = sign(a:: ) and 9_sign(a; ) = 25{x )): 

= 7^ / dy+\x+ - y+\j- (6) 



A^ 

in absence of background fields. Plugging Eq. ([6]) into Eq. ([4]) it turns out that there is a linear 
potential between charges. By dimensional analysis, g has dimensions of energy. 

We will consider this theory in the 't Hooft limit: large A and g held fixed. As A^ = 3 the 
color transformation property of a diquark (the scalar quark ipa) and of an antiquark are the same: 
diquarks are qq states with attraction in the color 3 channel. Extending the color group to SU {N) , 
with A^ > 3, this diquark-antiquark correspondence is lost. Here we treat the fields tfa as the fields 
of building block, pointlike, diquarks. As for color transformations, this is appropriate for A = 3, 
whereas for larger A^ it is an extrapolation. Alternatively one could consider the Corrigan-Ramond 
large A^ limit In the latter case quarks and 'larks' are introduced, transforming as the N and 
N(N — l)/2 representations of SU{N) respectively. Larks (. are represented by antisymmetric tensors 
iap — —£f3a coinciding with antiquarks if A' = 3. A theory of only larks is equivalent to QCD. In the 
CR large A^ limit, a baryon is represented by qq£ whereas there are no color singlets made up of three 
quarks in the standard large N. 

In the large A^ limit only planar diagrams are relevant and quark loops are suppressed. Hence gluon 
lines are impassable barriers as there are no gluon-gluon interactions and any gluon line crossing would 
violate planar ity. 
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Eq. ([2]) is derived from the Bethe-Salpeter equation expressed by the foUowing diagrammatic rela- 
tion: 





r p 




r ^ p 



where the disks on the external legs indicate the dressed propagators. The shaded blobs represent 
the matrix element of the time ordered product of two scalar quark fields between the vacuum |0) 
and the meson state \B) Q. We will denote by -0 the Fourier transform of this matrix element. The 
Bethe-Salpeter equation is the Dyson-Schwinger equation for a 4-point Green function G where a 
bound state \B) occurs as a mass pole (the residue at the pole giving the bound state wave function) 
in G: 



G 



\B)(B\ 



(7) 



The simplicity of the previous diagrammatic equation is a consequence of the large N limit, 
represent the Dyson-Schwinger series as: 



Let us 



o: 



In the large N limit, the only non- vanishing contribution from the amputated functions represented 
by the shaded boxes is the first one: 
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Once the gluon propagator is found to be (see Q): 



(8) 



where P indicates the Cauchy principal value, it is not difficult to compute the quark dressed propa- 
gators D{p). This can bee seen by observing that: 



(9) 



One obtains (cfr. [3|): 



Dip) = 



2p- 



(10) 
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where A is an infrared cutoff. To derive the latter expression it is convenient to observe that: 

^P+7r—-^^^—r-— = -*^sign(p_) (f 1) 



2p+p_ - A 

as can be proved by using the Plemelij identity 



/ 



dx = Idx-^^^^±nrfixo) (12) 

x-xoTie J (x - Xq) 



and that P/a; = 0. 

According to Eq. (llOp . in the hmit A ^ the quarks have infinite self-energy, hence are removed 
from the spectrum. Anyway the dependency on A is canceled in the Bethe-Salpeter equation which 
eventually defines the eigenvalue equation for scalar quark bound states; in other words we have a 
discrete spectrum of bound states but no free quarks. The diagrammatic equation given above writes 
as: 

i^{p, r) = g^D{p)D{p ^ r) j ^-^i^{P + k, r){2p + fc)_(2p - 2r + (13) 

Note the derivative couplings of scalar chromodynamics at the vertices of the gluon propagator. 
Observe also that the kernel of the principal value integral depends on minus momenta only: this is 
clear from Eq. ([6]) where only j_ currents are coupled to the linear potential, or equivalently from 
the expression of the gluon propagator. We can therefore integrate both sides of Eq. (fT^ over p+. A 
straightforward residue calculation gives: 

, 1 2^»e(p_)g(r_ -p_) 

''^+(p+-A)(p+-r+-i3)~ A-B^r+ ^'^> 

where we required a negative imaginary part for the pole A and a positive one for the pole B. If 
both poles were on the same half-plane the residue integration would yield zero closing the contour 
in the other half- plane. We close the contour of integration in the lower- half plane (clockwise sign). 
Defining: 

4'iP-,r)= / dp+ip{p+,p^,r) (15) 



/ 



one can write Eq. (|T3|) as: 

g^ %-)0(r_-p_) f- dfc_ ^ (p_+fc_)(p_+fc_-2r_) 

nP-^r)^-- — — / -^(/>fc_,r ^ (16 

27r r+-(A-B) (fc_-p_)^ 2p_ 2(p_-r_) 

The latter integral is infrared divergent when p^ ^ k^. We can expand it as: 

dfc_ f' dk^ , , dk_ 

(fcT^^ ^ Jo (A-^^"^(^-'^^l"J^oi__, (fc-:^ = 

This expression allows to cancel the 1/A infrared divergence in Eq. (fT6)) . Consider in fact that: 

TT A 2p_ 2(p-r)_ ^ ' 

Eq. ([2]) is obtained by defining r+ = ^"^ /{2r^), xr^ = p_ and t/r_ = fc_ and rescaling all the square 
masses in units of g^/vr. Observe that the ^-functions in Eq. ([T4|) define an interval outside which 
(f) = Q. This interval is p_ G [0, rJ\ or x G [0, 1]. In particular we have (/)(0) = </>(!) = 0. 

The integral in Eq. ^ gives its main contribution if y ~ x, where the kernels of Eq. ([T|) and ^ 
are the same. As the highest part of the spectrum is concerned, one can neglect Mf^s and the eigen- 
functions are approximated by 4'{x) ~ sinwa; = sinnTrx, with n integer n > 1. This approximation 
respects the condition (/)(0) = (t>{\) = 0. For periodic functions <j>{x) we have: 







exp{iujx) I exp(iwa;) . f exp(iwx) 
'^Vl ^ - 7 '^Vl ^ = f dy- = ~tt\uj\ exp(«wa;) (19) 
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as can be seen easily by applying p2|) . This result means that the eigenvalues in the highest part of 
the spectrum are /i^ ~ nn^, i.e., there is no continuum in the spectrum Q. The alternating parity 
for these 4>n{x) states is evident. 

In what follows we will focus on the lower part of the spectrum, and we will solve Eqs. ^ and ^ 
with a numerical approach. 

Discrete Spectra. Equations (|l|2p are integral equations with singular kernels and a prescription, 
the Cauchy principal value, on how to treat the singularity. In this section we will indicate the steps 
to put them in a form amenable to numerical computation. We follow a procedure which has been 
applied in Q (and reference therein) to solve Eq. ((!]). The main point is to write: 

N 

(f>{e) = amsinme (20) 

m— 

and set: 

1 + cos 6* 



(21) 
1 + cos 61' 

y = — ^ — (22) 

Integrating by parts in Eq. ([T]) to lower the singularity of the kernel, and exploiting the conditions 
(/)(0) = (/)(1) = 0, one obtains an integral of the form: 

,^osm^ (23) 
cos 6*' - cos 6* ^ ' 

which can be solved using the following relation between Chebyshev polynomials @ 

' dy- ^I'^y} ^ = 7rC/„_i(a;), (24) 



{y-x)y/l-y'^ 

where x S [—1,1]. Indeed this is done by observing that T„(cos6') = cosnO whereas J7„(cos6') 
sin(n + 1)6/ smO. The next step is to discretize the angle 6 according to the prescription jS)]: 



N + l 

and to use the following orthogonality relation: 

N N 



(25) 



sin k6„i sin kOn = ^ sin 9km sin 9kn = ^ t ^mn (26) 



fe=i fe=i 
to transform Eq. ([1]) into an eigenvalue equation of the form: 



N 

^ /n(^)„ _ ,,2 

m— 1 



O^^iam = f^f^^an (27) 



where we have 



,^ _ _ 4 sin 6lfc„ sin 6'fcm ( 2Af| 



A* + 1 smOk \ smOk i 



k=l 



The operator Onm can be diagonalized and a set of discrete eigenvalues /i(^) can be found. We aim 

here to find an operator oim and an eigenvalue equation form for Eq. ([2]) similar to that in ([27]) . The 
procedure to follow is basically the same as the one outlined above, just requiring some more algebra. 
An helpful relation to use to obtain Omn is provided by the following identity relating Chebyshev 
polynomials, similar to that given in Eq. ([M)) [l7| : 



' d,^iIZ^^ = -.T„(x) (29) 

-1 (y-x) 
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The result is: 



N r 



N - 



k=\ 



7 1 

TOTT — 

16 smt'fc 



sm ty„fc sm Ojnk 



— cos f^fc sm 6nk cos fc'mfe- 



-— cot6'fc sin6'„fc (sin6'mA: cos6lfc + sin 6'(„j+i)fc + sin 6'(„_i)fc) 

o 

mTT sin Onk 
32 sin6'fc 



sin( 



(m+2)fe 



sin( 



{m-2)k, 



(30) 



In the derivation of oim the trigonometric Werner identities have been used. Thus we have to pick 
up a minimum value of Mg to avoid tachyons. 

The operator obtained allows to write also Eq. ^ in a form amenable to numerical computation: 



N 



(31) 



and can be diagonalized to obtain the spectrum of eigenvalues A^(cr)- 

Numerical results. In Figure [1] we show a set of values (Ms, My) which give a ratio fi(a-)/ fJ,(Tr) 
compatible with the experimental one, taking the mass value of the a from py]. These values can be 



constituent quark model 



FIG. 1: The curve represents the {Ms,Mf) values which give the ratio /X(CT)/M{7r) of the lowest eigenvalues 
equal to the experimental one. We check that for large values of (M/,Ms) the two equations ([1]) and ((2]) are 
the same and the results converge with what expected from naive constituent quark model. The numerical 
analysis is performed by setting the A*' parameter in Eq. (|20|) to A'^ = 100. The resuhs are stable under 
variation of A'^. 




used to determine the gap between the ground levels and the first excited levels, having the desired 
parity, in scalar and spinor chromodynamics. The first radially excited 0^ state is expected to be the 
second excited state in the qq spectrum, the ground state being 0". Define the following ratio: 

(32) 

where with the superscripts 0,1,2,... we label the levels in the spectrum: is the ground state. 
The predicted value for the first positive parity radial excitation of a system of two scalar quarks (a 

[2] 

diquark-antidiquark state in our language) is then "i^J^ which can be obtained by: 

[2] m^' 1 

"^ J^ = T (33) 

using the mass m^/ — 1300 MeV for the first radial excitation of the pion and TOtt = 135 MeV, 
rria- = 470 MeV. With the values {Ms,Mf) shown in Fig. [Jit turns out that m[^j > 2590 MeV 
indicating that a radial excitation of the lowest lying tetraquark state would be likely a very broad 
object. Also there is no harm to find an exotic radial excitation in the range of masses where the 
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super-GeV scalar nonet is located. Similarly we can also estimate the lower bound for the first radial 
0~ excitation of a qq particle, finding m^jj^ > 3060 MeV. This can be done by assuming that the 
first qq scalar radial excitation is the /o(f370). The given minimum value occurs e.g. at the point 
{Mf, Alg) ^ (0.1, 4.3) of the parameter space spanned by the curve in Fig. [T] 

Heavy-light diquarks. Diquarks have been used in the literature also to discuss a number of 
recently discovered charmonium-like resonances, known as X, Y, Z particles, with phenomenological 
properties conflicting with those of standard charmonium states. The X(3872) found by Belle 
was recognized from the start to be a very anomalous charmonium state, despite its discovery decay 
mode into J/tpp that could have seemed the footprint of a higher charmonium. In JJJ the X{3872) is 
interpreted as a qq state with q — [cq], q being a light quark q — u,d. The diquark model predicts two 
charged almost degenerate X particles and two neutral ones, namely [cm][c'!2] and [c(i][cd]. As discussed 
in there is strong evidence by BaBar and Belle that a second neutral state, the X(3876), exists! As 
for now, there are no conclusive indications for X^. On the other hand Belle has found three charged 
resonances decaying into charmonium+charged tt, namely Zt (4051), .^2(4248), Z(4430). These are 
produced in B — > KZi, Z2 or Z. One can find a brief account on this in [W^ . 

For example, the ^(4430), observed in Z ^ V(25')7r+ [H, if confirmed, should necessarily be a 
multiquark object. In [T^ this state is interpreted as a [cu]i[cd]o radial excitation of a 1^ particle 
predicted in [llj and decaying into 'ip{lS)n^. 

Interestingly the difference in mass between the radially excited tetraquark Z+ and its fundamental 
state predicted in [11] is very close to M^(2S) — M^(i5). 

The notation [cw]i[cd]o means that we assume one of the two diquarks to have spin one. Light 
diquarks with spin one are thought to be less stable than spin zero ones, but, as a heavy-quark is 
introduced in the diquark, spin-spin chromomagnetic interactions are suppressed by the heavy mass 
and one expects to treat heavy-light diquarks with spin zero and spin one on the same footing. For a 
phenomenological investigation on higher light tetraquarks see [16]. 

In this section we want to investigate the system qiq -I- qqi using the methods developed above. 
We assume that the spin one diquark can be described by a colored (axial) field Ai in the 3-color 
representation. The gluon field in the adjoint is Ap with A] = and A* = —A*-' if the minimal 
coupling is ^ + gA^. We define J^^^ = [d^A], — d^A^^) for the axial diquark. 

The vertex AAA can be extracted by: 

-l=[AtA''^{F,,)) + A': TIM")) + T.M^^A'^))] (34) 
or in terms of Feynman rules: 



fc, a 

= -7^[ga^{k-p + r)^+ga^,{-p + r-2k)^+g^^{2p-2r + k)a] 



p — r + k,^ p — r,v 

where a,fi,iy = +, —. The light cone gauge fixes a = — and A = (e, — e). If the incoming ^ is ^ = (0, e) 
then the outgoing one is ^ = (e,0), i.e., at the vertices we have A+A- or vice- versa. We sum over 
these two alternatives. Taking a = = fi = +, from we get (— p+r — 2fc)_ + (2p — 2r + fc)_ = 
(p — r — k)- whereas taking a = +, /i = — we have {k — p + r)^ + (2p — 2r + fc)_ = {p~r + 2k)^. 
The sum of the two contributions involves (2p ~ 2r + k)^ as in (|13p . The same considerations are 
applied to the calculation of the self-energy (^^^) which therefore proceeds in the same way as in 
the scalar chromodynamics case. 

This allows to use the same Bethe-Salpeter equation given in ([2]) to study the case of charmed 
diquark- antidiquark states of the form [cu]i[c(i]o -I- [cu]o[c(i]i. 

In pT|], using chromomagnetic interaction Hamiltonians, the following mass values of 1^ , [eg] [eg] 
states, where found: Z(3754) and ^(3882). In |l3| the hypothesis is made that the observed 
Z+(4430) could be a radial excitation of the Z{3882). Here we assume that these three states, 
Z(3754), Z(3882), Z(4430) correspond to , /^I^l , of Eq. ©. If we require that /il^l/^M = 1.08 = 
Mz(3882}/Mz{3754) we See that the minimum allowed value of the parameter Ms (we could call it Ma 
in this case making reference to axial diquarks) is Ms ~ 10.43, quite higher than what found before. 
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consistently with the idea that we should describe spectra of charmed {i.e. heavier) particles. With 
the latter value of Ms one finds Ail'^V^'^' - 1-20 which requires that the value of the mass of Z(4430), 
as found in this model, is indeed 4658 MeV. This value is rather stable if one increases the parameter 
Ms up to Ms ^ 20. The next radial excitation is predicted at 6055 MeV. 

Conclusions. In this short note we assume that scalar chromodynamics in (l+l)-dimensions and 
with a large number of colors could be used to estimate the mass of the first radial excitation of a 
diquark-antidiquark meson. We applied a numerical method to solve the Bethe-Salpeter equations and 
compute the bound state discrete spectrum of this confining theory. The possible masses of the spinor 
and scalar quarks are found by imposing that the ratio of the ground state eigenvalues of the spinor 
and scalar Bethe-Salpeter equations, Eqs. ([T]) and ((2]) respectively, is equal to the ratio of the physical 
masses m^/mcr. Furthermore, with these masses, we were able to extract a minimum value for the 
first radial excitation of the ground state diquark-antidiquark spectrum. We extend our discussion 
to the heavy-light diquark sector finding that the Z+(4430) observed by Belle could correspond to 
the second radial excitation of the spectrum Z(3754), Z(3882), Z(4658)... where the lower states were 
predicted in 
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